Tidal effects around higher-dimensional black holes 
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In four-dimensional spacetime, moons around black holes generate low-amplitude tides, and the 
energy extracted from the hole's rotation is always smaller than the gravitational radiation lost 
to infinity. Thus, moons orbiting a black hole inspiral and eventually merge. However, it has 
been conjectured that in higher-dimensional spacetimes orbiting bodies generate much stronger 
tides, which backreact by tidally accelerating the body outwards. This effect, analogous to the 
tidal acceleration experienced by the Earth-Moon system, would determine the evolution of the 
binary. Here, we put this conjecture to the test, by studying matter coupled to a massless scalar 
field in orbit around a singly-spinning rotating black hole in higher dimensions. We show that 
in dimensions larger than five the energy extracted from the black hole through superradiance is 
larger than the energy carried out to infinity. Our numerical results are in excellent agreement with 
analytic approximations and lend strong support to the conjecture that tidal acceleration is the rule, 
rather than the exception, in higher dimensions. Superradiance dominates the energy budget and 
moons "outspiral" ; for some particular orbital frequency, the energy extracted at the horizon equals 
the energy emitted to infinity and "floating orbits" generically occur. We give an interpretation 
of this phenomenon in terms of the membrane paradigm and of tidal acceleration due to energy 
dissipation across the horizon. 
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I. INTRODUCTION 

Gravitational binaries are intrinsically complicated 
systems which display a wealth of interesting effects. One 
important effect, which occurs for example in several 
planet-moon systems, are tides generated by differential 
gravitational forces. Tidal forces in the Earth-Moon sys- 
tem have long ago locked the Moon in a synchronous 
rotation with the Earth and have increased the Earth- 
Moon distance P, [||. These processes, so-called tidal 
locking and tidal acceleration, respectively, are only pos- 
sible if there is some dissipation mechanism in the system. 
Because of friction, tides extract energy from the binary 
and, since angular momentum must be conserved, this 
provides a mechanism to exchange angular momentum 
between the Earth and the Moon. In the Earth-Moon 
case the dissipation is caused by the friction between the 
oceans and the Earth surface. In binaries containing ro- 
tating black holes (BHs) the event horizon can play the 
role of a dissipative membrane and tidal acceleration is 
known for many years under a different name: superra- 
diance 

Consider a Kerr BH whose angular velocity of the hori- 
zon is fin- A wave with frequency u> < mfln scattering 
off the BH (where m is the azimuthal quantum number) 
is amplified, and it extracts rotational energy from the 
BH. Superradiance is responsible for many interesting ef- 
fects (e.g. [ohTiI]). one of them is the possible existence 



* richard.brito@ist.utl.pt 
t vitor.cardoso@ist.utl.pt 
t paolo.pani@ist.utl.pt 



of "floating orbits" around BHs. Generically, orbiting 
bodies around BHs spiral inwards as a consequence of 
gravitational-wave emission. When the condition for su- 
perradiance is met, it is possible to imagine the existence 
of floating orbits, i.e., orbits in which the energy radi- 
ated to infinity by the body is entirely compensated by 
the energy extracted from the BH [5], 0|. Within gen- 
eral relativity in four dimensions, tidal effects are in gen- 
eral completely washed-out by gravitational-wave emis- 
sion and orbiting bodies always spiral inwards @. How- 
ever, when coupling to scalar fields is allowed, an induced 
dipolc moment produces a tidal acceleration (or polariza- 
tion acceleration [(|) which might be orders of magnitude 
stronger than tidal quadrupolar effects. Furthermore, in 
theories where massive scalar fields are present, the cou- 
pling of the scalar field to matter can produce resonances 
in the scalar energy flux, which can lead to floating orbits 
outside the innermost stable circular orbit [Tol . [l2| . 

It was recently argued via a tidal analysis framework 
that higher-dimensional BHs in general relativity should 
be prone to strong tidal effects @. One of the conse- 
quences of those studies was that orbiting bodies around 
higher-dimensional rotating BHs always spiral outwards, 
if the tidal acceleration (or, equivalently, the superradi- 
ance) condition is met. In this paper, we use a fully rela- 
tivistic analysis, albeit in the test-particle limit, to prove 
this behavior. For simplicity, we consider the coupling of 
massless scalar fields to matter around a rotating BH in 
higher-dimensional spacetimes. We show that, for space- 
time dimensions D > 5, tidal effects are so strong, that 
the energy extracted from the BH is greater than the en- 
ergy radiated to infinity. Higher-dimensional spacetimes 
arc of interest in a number of theories and scenarios El, 
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in our case we view them as a proof-of-principle for strong 
tidal effects in BH physics, without the need for reso- 
nances. We do not consider gravitational perturbations, 
gravitational effects should be subdominant with respect 
to the dipolar effects discussed here @. Nevertheless, the 
arguments presented in Ref. @ together with the present 
results show that a purely gravitational interaction also 
displays this phenomenon, which likely leads to new in- 
teresting effects in higher dimensional BH physics. 

This paper is organized as follows. In Section UH ex- 
tending the discussion of Ref. @ , we study the effect of 
an electrically charged particle orbiting a neutral central 
object in D = 4 + n dimensions. By applying the mem- 
brane paradigm [l4|, we derive a simple formula for the 
ratio between the energy flux at infinity and at the hori- 
zon in 4 + ro dimensions. In Section IIIII we compute the 
energy fluxes in terms of BH perturbations sourced by 
a test-particle in circular orbit around a spinning BH. 
We derive the Teukolsky equations and the expressions 
for the energy fluxes. In Section IIVI we solve the wave 
equation analytically in the low-frequency regime. We 
then compare our analytical results with those obtained 
by a direct numerical integration of the wave equation, 
as discussed in Section [V] We conclude in Section [VI] 
Throughout the paper we use G = c — 1 units, except in 
Sec. HIl where, for clarity, we show G and c explicitly. 



II. TIDES FOR CHARGED INTERACTIONS IN 
4 + n DIMENSIONS 

The flux emitted by a particle orbiting a spinning BH 
can be estimated at newtonian level in terms of BH tidal 
acceleration and by applying the membrane paradigm Q . 
In this section, we generalize the computation sketched 
in Ref. [6J to higher dimensions and to scalar fields. 

Let us consider the interaction of a particle with scalar 
charge q p and gravitational mass m p orbiting a neutral 
central object of mass M and radius R. If the object 
has a dielectric constant e = e, eo, the particle external 
field induces a polarization surface charge density on the 
central object and a dipole moment which are given re- 
spectively by [TH 



Cpol : 

P 



(3 + n)e l3E cos i? . 



(1) 

(2) 



charge q p due to the induced electric field is given by 

q P p 



sint9 . 



(4) 



Without dissipation, the dipole moment would be aligned 
with the particle's position vector. Here, we consider that 
dissipation introduces a small time lag r, such that the 
dipole moment leads the particle's position vector by a 
constant angle <f> given by (see @, for details) 



= (fi-n ff )r ; 



(5) 



where fin and fl are the rotational angular velocity and 
the orbital angular velocity, respectively. At first order 
in <j), the tangential component of the force reads 
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This exerts a torque r^F^ and the change in orbital en- 
ergy over one orbit reads 
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n(n - q h )t, 



(7) 



where, in the last step, we used Eq. ([2]). 

Remarkably, the equation above qualitatively describes 
the energy flux across the horizon of a rotating BH if one 
identifies iln with the angular velocity of the BH and 
the lag r with the light-crossing time, r ~ R/c, where 
R i+n _ G D M/[(n+l)c 2 }, G D is the D-dimensional grav- 
itational constant and M is the BH mass [iii . Accord- 
ingly, a particle orbiting a rotating BH in 4 + n dimen- 
sions dissipates energy at the event horizon at a rate of 
roughly 
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On the other hand, charged accelerating particles radiate 
to infinity according to Larmor's formula, which in 4 + n 
dimensions reads [161 ] (this will be also derived below, see 
Eq. @2)) 



where 



E 



Qp 



fyn+3)eoro + ™ 



(3) 



and ro is the orbital distance, Qr n +s) is the solid angle of 
the (ro + 3)-sphere, (3 is some constant that depends on 
the relative dielectric constant of the object and -d is the 
polar angle with respect to single axis of rotation of the 
central object. 

Assuming circular orbits, the tangential force on the 
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4+n 2 
r I 



(9) 



where 7 is some coupling constant. 

Tidal acceleration Q occurs when the orbit of the par- 
ticle is pushed outwards due to energy dissipation in the 
central object. This is only possible if two conditions are 
satisfied: (i) < Qh, so that Eh < and the energy 
flows out of the BH; and (ii) \Ejj\ > E OQl i.e. the rate 
at which energy is dissipated to infinity must be smaller 
than the rate at which energy is extracted from the BH. 
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From Eqs. ((5J and (|9]), and using il 
find 

\Eh\ _ K G ^ 



-(3+n)/2 



17f, we 



P (n + 3)(2-rQ 4 + r 



7^^nf^- n) -(c) ■ (10) 

where we have assumed £Ih S> Q and we have defined 
the orbital velocity 



v = [M(n + 1)] "+ 3 fl^ 



(11) 



At large distance, v ~ rofi and, when n = 0, we recover 
the standard definition, v = (Mfi) 1 / 3 . Surprisingly, for 
n > 1 (D > 5) tidal acceleration dominates at large dis- 
tances. This simple argument suggests that test-particlcs 
orbiting rotating BHs in dimensions greater than five 
would generically extract energy from the BH horizon 
at a larger rate than the energy emitted in gravitational 
waves to infinity. As a consequence, the orbital separa- 
tion will increase in time, i.e. the system will "outspiral" . 
In the next sections we shall prove this is indeed the case, 
by computing the linear response of a higher-dimensional 
spinning BH to a test-particle in circular orbit. 



III. SCALAR PERTURBATIONS OF 
SINGLY-SPINNING MYERS-PERRY BLACK 
HOLES 

A. The background metric 

In four dimensions, there is only one possible angu- 
lar momentum parameter for an axisymmctric space- 
time and rotating BH solutions are uniquely described by 
the Kerr family. In higher dimensions there are several 
choices of rotation axis, which correspond to a multitude 
of angular momentum parameters |l8j . Here we shall fo- 
cus on the simplest case, where there is only a single axis 
of rotation. In the following we shall adopt the notation 
used in Refs. (l9l - [2l| . to which we refer for details. 

The metric of a 4 + n dimensional Kcrr-Myers-Pcrry 
BH with only one nonzero angular moment um p arameter 
is given in Boyer-Lindquist coordinates by [18( 



A-a 2 sin 2 tf ji2 2a{r 2 + a 2 - A) sin 2 d 

S 

sin 2 tfdcf) 2 + —dr 2 



ds 2 = dt 2 

(r 2 +a 2 ) 2 - Aa 2 sin 2 tf . 



Y,dd 2 + r 2 cos 2 Mtil 



where 

£ = r 2 



J 2 cos 2 ■& . 



A = r 2 + a 2 - 2Mr x 



dtdcj) 



(12) 



(13) 



and dil 2 t denotes the standard line element of the unit 
n-sphcre. This metric describes a rotating BH in asymp- 
totically flat, vacuum spacetime, whose physical mass M. 



and angular momentum J (transverse to the r<\> plane) 
respectively read 



M= in + 2)A n+ , M 



n + 2 



Ma . 



(14) 



where A n+2 = 27r n+3 / 2 /T[(n + 3)/2]. 

The event horizon is located at r = rjj, defined as 
the largest real root of A. In four dimensions, an event 
horizon exists only for a < M. In five dimensions, an 
event horizon exists only for a < \/2M, and the BH area 
shrinks to zero in the extremal limit a — > \2M . On the 
other hand, when D > 5, there is no upper bound on the 
BH spin and a horizon exists for any a. 



B. Setup 

We consider a small object orbiting a spinning BH and 
a massless scalar field coupled to matter. At first order 
in perturbation theory, the scalar field equation in the 
background <jT2j) reads 



Dip = 



1 







-g cto 1 



(15) 



where a is some coupling constant. For simplicity we 
focus on source terms of the form 



T 



dr 
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m p S i4+n) (x-X(T)), 



(16) 



which corresponds to the trace of the stress-energy ten- 
sor of a point particle with mass m p . We also restrict to 
equatorial circular orbits (z? = 0, d — tt/2), which is an 
unrealistic approximation in higher dimensions: generic 
circular orbits are unstable, with an instability timescale 
of order of the orbital period [17j . Nevertheless, our pur- 
pose here is to show that tidal effects can dominate, it is 
not clear what the overall combined effect of tidal accel- 
eration and circular geodesic motion instability is. Ex- 
tending the present analysis to generic orbits and relax- 
ing the test-particle approximation are interesting future 
developments. 

For prograde orbits around a singly-spinning Myers- 
Perry BH (|12[) the energy, angular momentum and fre- 



quency of the point particle with mass m p orbiting at 
r = r read (l7| 



E,, 



a^(n + l)M + r 2 



2Mr' n 



2ay/(n + 1)M + r 2 - (n + 3)Mr : 
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y/(n+l)M (r 2 



2a-s/(n + 1)M + r 2 - (n + 3)Mr ; 



(18) 
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V(n + 1)M 
iy/(n+l)M + 



(19) 



The only nonvanishing components of the (4+n)-velocity 
XJ V of the particle on a timelike geodesic are given by 

f ( o 2A/a 2 \ 2MaL„ , . 

' ^"-^(20) 
'o / 'o 



m p A r = ro U l = lr z +a 2 + E. 



m p A r=ro l^ 



2MaE p 



2M 



l 



(21) 



3r 2 A 2 A [3r 2 + a 2 - 2A/r 1 -"(2 - n)] 

V= cj 2 + 



(r 2 + a 2 ) 



2^4 



- a 2 ) 2 { 



(r 2 + a 2 ) 
■ A 



a 2 m 2 - 4M + 



( r 2 +fl 2)3 

" - A(^ 2 a 2 + A 



,71— 1 



n(2-n)A 2n(l-n)M 

- n 4 



4r 2 

j(j + n- l)a 2 



2 r n+l 



C. The wave equation 

Because of the coupling to matter, the orbiting object 
emits scalar radiation which is governed by Eq. (|15|) . To 
separate Eq. (|15p . we consider the ansatz 

ip{t,r,d,4>)=Y. j ' due im *- iut R(r)S lm M Y ^ ( 22 ) 

l,m 7 j 

where Yj are hyperspherical harmonics [2(| [HJ on the 
n-sphere with eigenvalues given by —j(j + n — 1) and 
j being a non-negative integer. The radial and angular 
equations read 



dr \ dr 
j{j + n- l)a 2 



[w(r 2 + a 2 ) — ma] 



-\}R = T h 



(23) 



and 



d f . _ a dSl rn j 

sin $ cos v - 

dv 

j(j + n- 1) 



sin $ cos" $ di9 



[w 2 a 2 cos 2 



sin 2 



cos 2 19 



+ 4 



iroj 



Slmj — . 



(24) 



where A = Ai rn j — 2muja + cj 2 a 2 and 

^=-^^•(^2)^(^/2^/2,...) 
x <5(r — 7'o)(5(mft p — cj) , 



In the low frequency limit the angular equation 
can be solved exactly. At hrst order in au, the eigenvalues 
can be computed analytically [22j 

A kjm = {2k+j + \m\){2k+j + \m\+n+l) + 0{auj) . (28) 

By setting 2k = I — (j + \m\), the eigenvalues above take 
the form Aij m = 1(1 + n + 1) and I is such that I > 
(j + \m\), which generalizes the four-dimensional case. 
An important difference from the four-dimensional case 
is that regularity of the angular eigenfunctions requires k 
to be a non-negative integer, i.e. for given j and m only 
specific values of I are admissible. In fact, it is convenient 
to label the eigenfunctions and the eigenvalues with the 
"quantum numbers" (k, j, m) rather than with (l,j, m) as 
in the four dimensional case. The (non-normalized) Oth- 
order eigenfunctions are given in terms of hypergcomctric 
functions [H H3 

■i , n+1 . n+1 , 
-k,k+j + \m\ + — — — ;x* 



S kjm oc sin(tf)l"V> 



(25) 



where x = cos(#). We adopt the following normalization 
condition 

r-7r/2 

dd Sin ■& COS™ tiSkjm Skjm = 1 » (29) 

where the integration domain has been chosen in order 
to have a nonvanishing measure also in the case of odd 
dimensions. Note that this normalization differs from 
that adopted in Ref. $2%. 

We note that at ■& = tt/2 only hyperspherical harmon- 
ics with j = are non- vanishing. Thus, in order to cal- 
culate the fluxes on circular orbits, one only needs to 
consider terms with j = 0. In this case, the hyperspheri- 
cal harmonics Yq are constant. 



which has been derived from the stress-energy tensor 
of the point particle. Defining a new radial function 

x lmj{r) 

X lmj =r n / 2 (r 2 +a 2 y/ 2 R : (26) 
we get the non-homogeneous equation for the scalar field 



dr 2 



V 



X lmj (r*) 



A 



(r 2 + a 2 )3/2 



„n/2j. 



Imj 



(27) 



where dr/dr* = A/ (r 2 + a 2 ) defines the standard tortoise 
coordinates and the effective potential V reads 



D. Green function approach and energy fluxes 

To solve the wave equation, let us choose two inde- 
pendent solutions X r k " m and X^ m of the homogeneous 
equation which satisfy the following boundary conditions 



kjm 



kjm ■"■out^ i -fiinC 



as r — > oo . 



yoo 
kjm 



B out e lkHr ' + B- m e- ikHr * 



X 



kjm 



as r — > Tu ■ 



(30) 
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Here ku = u> — mfijy, fcoc = w and 17// = 
-lim r ^ rfJ g t< t,/g^ = a/(r^ + a 2 ) is the angular ve- 
locity at the horizon of locally nonrotating observers. 
The Wronskian of the two linearly independent solutions 
reads 



W = xri 



dX 



kjm dr* 



-X 



klm dr, 



2ikooAi n , (31) 



and it is constant by virtue of the field equations. Finally, 
Eq. (|27j) can be solved in terms of the Green function |23j 



Xkjm [t* ) 

X th 

L kjm 



^kjm 

w 



X T , H dr 



+ 



Ar'™/ 2 



Tkjm(r) , /2 + a2 ^ 3 / 2 X k] 



(32) 



For very large values of r the above equation has the 
following asymptotic form 



X k j m (r oo) 



T kjm (r')X% 



Ar'™/ 2 
— drL 



( r /2 + a 2)3/2 



= Z% m S{u> - mfy e ifc ^« , 
where, using Eq. (|25() . 

^fm( r o)^. m (7r/2)y;(7r/2,7r/2,...) 



J kjm 



V r o+ a2r o 



n/2 



(33) 



(34) 



Likewise, at the horizon we get 

X kjm {r* -»■ -oo) = Z£f m tf(w - mf2 p )e- ife ^ r * (35) 
where, 



kjm 



Xkjmiro) S* kjm (ir/2)Y?(n/2, n/2, . . .) 
-a — ; — 7- m r 



WW 



V r o + a2r o 



n/2 



(36) 

The scalar energy flux at the horizon and at infinity 
arc defined as 



n 

Eh. oo = hm / ddd<f> I I d$iy/—gT t 

r—^r/f ,oo / J- 

2—1 



where the stress tensor reads 



= (V^V^* - -g ul ,V a ipX7 a ip* 
Finally, using Eqs. and ([3"5"j) . we get 



fejm 



rH-,oo |2 



(37) 



(38) 



(39) 



The equation above shows that, if the superradiant con- 
dition kn < (uj < m^ln) is met, the energy flux at the 
horizon can be negative, Eh < 0, i.e. energy can be ex- 
tracted from a spinning BH [f| Hj]. In four dimensions, 
\Eh\ -C Eoo and the superradiant extraction is generi- 
cally negligible. As we show in the next section, in higher 
dimensions the opposite is true, \Ejj\ S> E^ and super- 
radiance dominates over gravitational-wave emission. 



IV. ANALYTICAL SOLUTION AT LOW 
FREQUENCIES 

The scalar flux can be evaluated analytically in the 
low- frequency regime (see e.g. (2514291 ] ) . Let us first focus 
on the solution X T ^ m , which is regular at the horizon. 

We first make the following change of variable 



h = 



A 



r- + a 



dh , , . A(r) 



dr 



r* + a z 



(40) 



where A(r) = (n + 1) + (n — l)a 2 /r 2 . Then, near the 
horizon r ~ rjj, the radial equation (|23|) can be written 

as 

, . , . d 2 R , — j . dR 
h{1 - h) ^ + {1 - D * h) dk + 



P 2 



A 



R = 0, (41) 



A(r H ) 2 h{\ - h) r 2 H A(r H ) 2 (l-h) 
where 

P = oj(r H + a 2 /r H ) - ma/r H , 
A = [1(1 + n + 1) + j(j + n - l)a 2 /r 2 H ](r 2 H + a 2 ) , (43) 



(42) 



1 



^a 2 r 2 H 



[{n+l)r 2 H + {n- l)w 



n 2 • 



(44) 



Using the redefinition R(h) = h a (l - h) t3 F(h) the above 
equation takes the form 



d 2 F dF 
k{1 ~ h) lh 2 - + [c - {a + b+m ih 

with 



(ab)F = 0, (45) 



a = a + /3 + D*-l, b = a + /3, c=l + 2a, (46) 

and a and /3 must satisfy the following algebraic equa- 
tions 



P 1 



A(r H ) 



0. 



P 2 +P(D,- 2) + 

whose solutions read 
P 

a± = a 



p2 



A 



A(r H ) 2 r 2 H A{r H f 



(47) 



0, (48) 



At 



A(r H ) ' 

(2-0.: 



(49) 



±W(Z>*-2)= 



p 2 



A 



A(r ff ) 2 r H A(r H ) 



.(50) 



The two linearly independent solutions of Eq. (|45l) are 
i^(a, &; c; /i) and /i 1_c F(a+ 1 — c, 6+ 1 — c; 2 — c; /i), where 
F is the hypcrgcomctric function. Convergence requires 
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Rc [c — a — b] > 0, which can be only obtained if the mi- 
nus sign is chosen in the solutions above. In the follow- 
ing, we shall identify /3 = /3_ and a = a_. The general 
solution of Eq. ([4"Tj) is then 

R(h) = A^il - hfF(a, b, c; h) 

+ Bi/i" a (l - hfF{a + 1 + c, b + 1 - c, 2 - c; h). (51) 
Expanding the above result near the horizon, we get 

R(h) = Aih^^hn +B x ti^^ = A 1 e~ ikHr '+B 1 e lkHr *. 

(52) 

Regularity at the horizon requires B\ = 0. The near- 
horizon solution can be written as (see f30l ]) 

R(k)=A lh *(l-hf r[l + 2a]r[2-I>.-20 v 



r[2-D» + a-j8]r[l + a-jS] 



F(a, 6, a + 6 - c; 1 - h) + A x h a (l - h) 
T[l + 2u]T[2p + D* -2] 

rlQ + zS + D, - i]i> + /?] x 
F(c — a, c — b, c — a — b + 1; 1 — h). 



(53) 



We can now expand this result in the low-frequency 
regime and for small values of a/rji, in the region where 
1 — h <C 1 and r ^s> Th , 



R 



V r H 

J, 



T[l + 2a]r[2 - £>» - 2/3] 



(2M) 



2(™+i) rfv 1 r - iy * 



+ a - £]r[l + a - 0\ 



(54) 



We will follow Poisson to solve the wave equation at 
large distances [25|. It is useful to rewrite the radial 
equation (|27p in terms of the dimensionless variable z = 
ur. At large distances Eq. (|27p reads 



d 2 (n + l)e d 
+ 

[! + ■ 



f dz 2 



1(1 



l) + f(l + f) 



z 2 +™ 



1)] 



r 3+n 



^ 3 m(z) = , 



(55) 



where / = 1 — e/z" +1 , e = 2Mu n+1 is a dimensionless 
parameter, and we used the fact that at large distances 
the eigenvalues take the form Aij m = 1(1 + n + 1). We 
can rewrite it in a simpler form if we define the quantum 
number J(J + 1) = 1(1 + n + 1) + |(1 + §). Solving for 
J, and assuming J is a non- negative number, we get 



J = l + 



In the limit Eq. (|55j) reads 



| 1 JQ7 + 1) 
<iz 2 z 2 



X kjm (z) = . 



(56) 



(57) 



The solution can be written in terms of a linear com- 
bination of Riccati-Bcsscl functions, \fzJj+\/2( z ) arL d 



y/~zNj +1 / 2 (z). The requirement that X k " m (z) be regu- 
lar at the horizon demands 



X r k « m (z) = B^J, J+1/2 (z). 



(58) 



where B is a constant. The asymptotic expansions for 
the Bcsscl functions are well known and read 



Bz J+1 



2 J + 1 /2r[J + 3/2] 



[l + 0(z 2 )]. (59) 



At large distance, the second term within the square 
brackets is subdominant and we shall ignore it. Match- 
ing (HU) to ([HI) we get 



_ 2 J+1 ' 2 T[J + 3/2]r[l + 2a]r[2 - D* - 2/3] 
~ r[2 - £>* + a - 0\T[1 + a - 0\ X 

( 2M )l/(2n+2) 



E (J+l)/(n+l) r V 2 



[I + 0(e)] . 



(60) 



The parameter Aj n can be extracted from the behavior of 
the function near z = oo. Recalling the large- argument 
of the Bessel functions, X k " m (z — > oo) ~ By / 2/nsm(z — 
Jir/2) and, using Eq. (|30|) . it follows that 



A 5 , 



2 J r[J + 3/2]r[l + 2a]r[2 - D* - 2/3} 
0Fr[2 - D* + a - /3]r[l + a - /?] 



,l/(n+l) 



(2M) 



l/(2n+2) 



1/2 
"if 



[l + 0(e)]. (61) 



With all of this at hand, we can now compute the flux 
at infinity in the low-frequency regime. From Eqs. (|34|) 
and (|39"T) we get 



m 



m 2 n 2 p \z^ m \ 



2+2Z+n 



_2'+ n / 2 + 1 rp + n/2 + 3/2] 
[(n + l)M] l+n/2+1 \S kjm (7r/2)\ 2 \Y 3 (7r/2 7 n/2, . . .)| 2 



2i(n + l) + (,i + 2)(,i + 3) 



x r n 



(62) 



where we used the fact that for small frequencies (large 



distances) 



a ~ r 

-(3+n)/2 



u 1 



1 and uj 



m^J(n + l)Mr 

Let us now perform the same calculation for the solu- 
tion X^ m , which satisfies outgoing- wave boundary con- 
ditions at infinity. The method is analogous to that al- 
ready described above. However, since for this case the 
boundary condition is imposed at infinity we do not re- 
quire regularity at the horizon. In the limit e< 1, X^ m 
can be identified, up to a normalization constant, with 



Xkj m { z ) — C y/z Hj^/^z) , 



(63) 



where Hj^ x , 2 (z) — Jj+i/2( z ) + iN J+1 / 2 (z) is the Hankel 
function. To determine the constant C wc match this 
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solution in the limit z — > oo to the required boundary 
condition 



X£ m (z -> co) 



(64) 



Recalling the asymptotic behavior of the Hankel func- 
tions, we have 



CV^hW 1/2 (z oo) ~ CyJ ^e iz [(-i) J+1 + 0(l/z)} , 

(65) 

from which we get C = i J+1 y/n/2. 

The small-argument behavior of the Bessel functions 
reads (when J > 0) 



N J+1 / 2 (z « 1) ~ 



, .7+1/2 



2 J+i/2 r [j + 3 /2] 

2J+1/2 r[j+ 1/2] 



l + 0(z 2 )] 
L l + 0(z 2 )] 



(66) 



At leading order and near z = 0, the function N J+1 / 2 
dominates over Jj+i/%. Hence, we get 



^Jy TO (*<l)~i J -s r[J + l/2K 



(67) 



We can now compute the flux across the horizon. Using 
Eqs. (1311) and flM) we get 



n + 1 



1/2 



r H (2M)- 



1+1 X 



(68) 



|^ m ( 7 r/2)| 2 |y j -( 7 r/2 j7 r/2 ) ...)| 2 xr 2 

wllprp F _ r[l+n/2+l/2]r[2-D,+ a -f)]r[l+ a -P] 
wiicic x 1 2r[Z+ri/2+3/2]r[l+2a]r[2-_D»-2/3] " 

We can now obtain an expression for the ratio of the 
fluxes on the horizon and at infinity for general I, m 
and n. Using the expressions for Eh and E^ calculated 
above, we find 

E H _ k H r H [2'+»/2+i r [; + n/2 + 3/2]] 2 



,21+n+l 



n+1 



(2M)" 



(n-l)(2i + n+l) (n-l)(n+l+31) 







(69) 



This can be written as a function of the orbital velocity, 

I 2 



E H _ k H r H [2 / +"/ 2 + 1 r[Z + n/2 + 3/2]]' 



n + 1 



■nm 

1 + 21 + 



2l+n+l 

(n- l)(n+l+2Q 



rf x 



X V 



l+l 



(70) 



where we have used Eq. (jTTJ) at large distance. For suffi- 
ciently small orbital frequencies, such that the superradi- 
ance condition is met, the flux at the horizon is negative, 



we then find that the ratio between the fluxes grows in 
magnitude with ro and the particle is tidally accelerated 
outwards. For the dipolar mode, 1 = 1, this expression 
is in complete agreement with the expected behavior de- 
rived from a newtonian tidal analysis, cf. Eq. (|10[) . 

Note that these results were derived under the assump- 
tion of slow rotation, a -C th- This approximation is 
particularly severe in the near-extremal, five-dimensional 
case, where r# — > 0. Nevertheless, as we discuss in the 
next section, our method captures the correct scaling of 
the energy fluxes for any spin and it even gives overall 
coefficients which are in very good agreement with the 
numerical ones in the slowly-rotating case. This is shown 
in Fig. [TJ where we compare the analytical results of this 
section with the numerical fluxes computed in Sec. |V] 



10 



10 



10' 



10 l 



10' 



10 



D=4 
D=5 
D=6 
D=7 
D=8 




0.001 



FIG. 1. Comparison between the flux ratio p = E^ 1 /E™ 
(in absolute value) calculated analytically and numerically, 
as a function of the particle velocity v for n = 0,1,2,3,4 
(D = 4,5,6,7,8) and a = 0.2M 1/(1+n) . The straight curves 
correspond to the analytical formula with / = 1 and the dots 
are the numerical results discussed in Sec. [V] In the slowly- 
rotating regime, numerical results are in very good agreement 
with the analytical formula. 



V. NUMERICAL RESULTS 

The Green function approach described above can be 
implemented numerically using standard methods flol . 
Ell H3[. For given values of ro, a and n. we can com- 
pute the fluxes by truncating the sum in Eq. (|39[) to some 
^max; wimax and j max . As discussed before, for circular 
orbits only j = terms give a nonvanishing contribution. 

For small and moderately large orbital velocities, the 
sum converges rapidly even for small truncation orders 
and we typically set fc max = 3 and m max = 6. However, 
the convergence is very poor when the orbital velocity 
approaches the speed of light, i.e. when the orbit is close 
to the prograde null circular geodesic. Recall that circu- 
lar orbits around Myers-Perry BHs in higher dimensions 
arc unstable [l7| and, in particular there is no innermost 
stable circular orbit for D > 4. Thus, for our purposes we 
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could in principle consider particles in circular orbit up to 
the light- ring, which exists for any dimension [TtJ • As the 
particle approaches the light-ring, the flux is dominated 
by increasingly higher multipoles, thus affecting the con- 
vergence properties of the series (|39|) . For this reason, 
the plots presented below are extended up to some value 
of the velocity that guarantees good convergence. 

Furthermore, for large orbital velocity and highly spin- 
ning BHs, the zeroth order angular eigenfunctions and 
the corresponding eigenfrequencies ([2"5)l might be not ac- 
curate. Therefore, when aw > 1, we have used exact nu- 
merical values of Akj m obtained by solving Eq. (|24p with 
the continued fraction method |22[. We note however 
that Eq. (|28|) reproduces the exact results surprisingly 
well, even when au ~ 1. 

We checked our method by reproducing the results of 
Ref. [12j j for the massless case in four dimensions. In 
addition, we can compute the energy flux in any number 
of dimensions. The fluxes Eh and Eoo for D = 5 and 
D = 6 are shown in Tables U and |TT] for r = 10ry. 
We show the total flux as well as the first multipolar 
contributions. 

TABLE I. Fluxes across the horizon and to infinity for n = 1 
(D = 5), a = M 1/{1+n) and r = 10r H - In the last row we 
show the total flux obtained summing up to fc max = 3 and 

"imax — 6. 



k 


III 


J 


ro/r H 


E H (am p ) 2 




2 


[EhI/E^ 





1 





10 


-1.4759 x 10" 


-9 


1.5790 x 10 


-!) 


0.9347 





2 





10 


-8.3175 x 10" 


11 


1.6288 x 10" 


10 


0.5107 





3 





10 


-3.1691 x 10" 


12 


1.0080 x 10" 


-11 


0.3144 


1 


1 





10 


-1.2718 x 10" 


14 


5.1192 x 10" 


-16 


24.844 


^-jkmj 


10 


-3.1248 x 10" 


-9 


3.5052 x 10 


-9 


0.8915 



TABLE II. Same as in Table U but for n = 2 (-0 = 6). 



k 


m 


J 


ro/rn 


E H {ctm p )~ 2 


fioo(am p ) 2 


\Eh\IEoo 





1 





10 


-4.1768 x 10" 12 


1.3915 x 10" 14 


300.155 





2 





10 


-3.0159 x 10" 13 


3.8070 x 10" 16 


792.180 





3 





10 


-1.3677 x 10" 14 


4.8651 x 10" 18 


2811.35 


1 


1 





10 


-2.8126 x 10" 16 


1.0497 x 10" 22 


2.6796 x 10 6 


2— /fcm7 


10 


-8.9858 x 10" 12 


2.8601 x 10" 14 


314.176 



TablcsHlandHIlconfirm our analytical expectations that 
the behavior for n > 1 (D > 5) is qualitatively different: 
the energy flux across the horizon is larger (in modulus) 
than the flux at infinity. This is shown in Fig. [21 where 
we compare the flux ratio p = E H ot /E™ as a function 
of the orbital velocity v for a = 0.99M 1 /( 1 +") in vari- 
ous dimensions. This figure is analogous to Fig. [T] but 
for a = 0.99M 1 /( 1+n \ i.e. a regime that is not well de- 
scribed by the analytical formula (|T0|) . For D = 4, we 
find the usual behavior, i.e. the flux at the horizon is 
usually negligible with respect to that at infinity and the 
ratio decreases rapidly at large distance. The case D = 5 



marks a transition, because p is constant at large dis- 
tance. This is better shown in the left panel of Fig. [3] 
On the other hand, for any D > 5 the flux across the 
horizon generically dominates over the flux at infinity. 



-» 1 1 








D=4 






-- D=5 






■— D=6 






D=7 






D=8 














-I 


— | £1=0.99 M 1,(1 * n) 





0.001 0.01 0.1 



V 

FIG. 2. The flux ratio p = E"£*/E™ (in absolute value) as 
a function of the particle velocity v defined in Eq. (|lll) for 
n = 0, 1, 2, 3, 4 (D = 4, 5, 6, 7, 8) and a = 0.99Af 1/(1+ "' . 

In Fig. [3] we show the flux ratio p for some selected 
value of the spin parameter a in five dimensions (left 
panel) and in six dimensions (right panel). When D = 5, 
the ratio is constant in the small v region and it ap- 
proaches unity in the extremal limit, a — > s/2M. As 
shown in the right panel of Fig. [31 when D = 6 there ex- 
ist some orbital velocity for which —p = l, corresponding 
to a total vanishing flux, Eh + E^ = 0. These orbital 
frequencies correspond to "floating" orbits 0, [To[ ■ Al- 
though in the right panel of Fig. [3J this is shown only 
for a/M 1 / 3 = 0.1, 0.2, 0.3, we expect this to be a generic 
feature also for larger values of the spin. The poor con- 
vergence properties of the series (|39[) prevent us to extend 
the curves to larger values of v, where floating orbits for 
a > 0.3A/ 1 / 3 are expected to occur. At smaller velocity, 
the energy flux contribution dominates and the motion of 
the test-particle is generically dominated by tidal accel- 
eration. Similar results can be obtained for any D > 6. 



VI. CONCLUSIONS 

We computed the rate at which the energy is extracted 
from a singly-spinning, higher-dimensional BH when a 
massless scalar field is coupled to a test-particle in circu- 
lar orbit. We showed that, for dimensions greater than 
five and small orbital velocities, the energy flux radiated 
to infinity becomes negligible compared to the energy ex- 
tracted from the BH via superradiance. 

Although we considered scalar-wave emission, we ex- 
pect our results to be generic in higher dimensions. In 
particular, superradiance should be a dominant effect 
also for gravitational radiation. At leading order, the 
ratio \Ejj\/Eoo for gravitational radiation should scale 
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FIG. 3. The ratio p = E^/E™ as a function of the orbital velocity defined in Eq. for several values of a. Left panel: 
when D = 5, the ratio is constant in the small v region and it approaches unity in the extremal limit, a — > \/2M. Right panel: 
when D — 6, the flux at the horizon can exceed the flux at infinity. For each curve, the intersection with the horizontal line 
corresponds to a floating orbit, —p = 1. Note that, at large orbital velocity, the superradiant condition is not met and Eh > 0. 



with the velocity as described by Eq. (|70p . The domi- 
nant quadrupole term (I = 2) reads [6| 

\Eh\ («-l)(«+B) 

'— L V " + 1 

Erx, 

By comparing the formula above to Eq. (|70p with 1 = 1, 
we note that dipolar effects are dominant over their 
quadrupolar counterpart. Nevertheless, even in the 
purely gravitational case, tidal acceleration and floating 
orbits around spinning BHs are generic and distinctive 
effects of higher dimensions. 

In principle, gravitational waveforms would carry a 
clear signature of floating orbits [HI [H| . Does floating or 
these strong tidal effects have any significance in higher- 
dimensional BH physics? We should start by stressing 
that circular geodesies in higher dimensions are unstable, 
on a timescale comparable to the one discussed here [ItJ ; 
however, our analysis suggests that, while more pro- 
nounced for circular orbits, tidal acceleration is generic 
and in no way dependent on the stability of the orbit 
under consideration. We are thus led to conjecture that 
tidal effects are crucial to determine binary evolution in 



higher dimensions. It is possible that tidal effects already 
play a role in the numerical simulations of the kind re- 
cently reported in Rcfs. [l3|, HH, [32|, but further study 
is necessary. One of the consequences of our results for 
those type of simulations is, for instance, that in higher 
dimensional BH collisions the amount of gravitational ra- 
diation accretion might play an important role. It would 
certainly be an interesting topic for further study to un- 
derstand tidal effects for generic orbits, and to include 
finite-size effects in the calculations. 
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